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ABSTRACT 

Observations of extragalactic objects need to be corrected for Galactic absorption and this is often 
accomplished by using the measured 21 cm HI column. However, within the beam of the radio telescope 
there are variations in the HI column that can have important effects in interpreting absorption line 
studies and X-ray spectra at the softest energies. We examine the HI and DIRBE/IRAS data for lines 
of sight out of the Galaxy, which show evidence for HI variations in of up to a factor of three in 1° fields. 
Column density enhancements would preferentially absorb soft X-rays in spatially extended objects and 
we find evidence for this effect in the RO S AT PSPC ohsevv&iions of two bright clusters of galaxies, Abell 
119 and AbeU 2142. 

For clusters of galaxies, the failure to include column density fluctuations will lead to systematically 
incorrect fits to the X-ray data in the sense that there will appear to be a very soft X-ray excess. This 
may be one cause of the soft X-ray excess in clusters, since the magnitude of the effect is comparable to 
the observed values. 



Subject headings: ISM: structure - X-rays: galaxies: clusters: ISM 



1. INTRODUCTION 

In the course of analyzing X-ray observations of extra- 
galactic objects, corrections are applied for the effects of 
Galactic absorption. The amount of Galactic absorption 
can be fit directly from the observations for adequately 
strong sources with a known underlying spectrum, and this 
value can be compared to the amount of HI measured from 
' 21 cm line emission. The similarity or difference between 
these two measures of the absorption column can be used 
for a variety of purposes, such as to examine whether there 
is excess absorption within a cluster of galaxies (White 
et al. 1991; Arabadjis & Bregman 2000) or whether ab- 
sorption by molecular material is commonplace within the 
Galaxy (Arabadjis & Bregman 1999a). In other cases, one 
can fix the X-ray absorption column at the 21 cm HI value 
and examine whether there is an additional emission com- 
ponent in a system, such as the very soft X-ray excess (the 
EUVE excess) that has been claimed to exist in clusters 
of galaxies (e.g.. Lieu, Bonamente & Mittaz (1999); but 
see Bowyer, Berghofer & Korpela (1999) and Arabadjis & 
Bregman (1999b)). 

In these studies, one usually ignores the small-scale 
structure of the neutral Galactic ISM, assuming that the 
mean 21 cm column density within a radio telescope beam 
is uniform across the face of the beam, which is typically 
30' for the single-dish radio telescopes that most of the HI 
values are based upon. Naturally, there is structure to the 
total neutral column within the radio telescope beam and 
this has the potential of infiuencing the analysis of extra- 
galactic X-ray emission as well as absorption studies in the 
optical and ultraviolet region. Here we examine the mag- 
nitude of the variation of the small-scale structure in the 
neutral Galactic layer and we discuss the implications for 
the analysis of data. First, we present an idealized model 



to illustrate the magnitude of the effect (§2), and then 
we present the HI and X-ray observations that show the 
magnitude of the variation in the absorption column (§3, 
4). Among other effects, we show that by ignoring such 
structure in the ISM, absorption corrected X-ray spectra 
of clusters will systematically show an apparent excess soft 
emission component, of the sort that is seen in some EUVE 
observations of galaxy clusters (§5). 

2. AN EXAMPLE OF THE IMPORTANCE OF SMALL-SCALE 
STRUCTURE 

The importance of small-scale structure on absorption 
corrections can be illustrated by adopting simple models 
for the fiuctuations in Nh that are sensible and are consis- 
tent with the observations discussed below. In one model, 
we assume that the column density distribution f{N) has 
a minimum level, Niow, a maximum level, N^igh, and be- 
haves as a power-law between these levels so that 
/(A) = A-™ for Nio^ <N< N^gh 
f{N) =0 for A < Aio„, Nh^gh < N (1) 

The mean column density is the 21 cm HI column, which 
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and in the case of Nhigu — oo, we have the simple form 
A2icm = ^E^Niow For individual lines of sight out of the 
Galaxy, one is more likely to measure a column that is typ- 
ically closer to the median column density rather than the 

^ University of Michigan, Department of Astronomy, 833 Dennison Building, Ann Arbor, MI 48109 
^ Institute for Astronomy, University of Hawaii, 2680 Woodlawn Drive, Honolulu, HI 96822 

^ Massachusetts Institute of Technology, Center for Space Research, 70 Vassar Street, Room 37-501, Cambridge, MA 02139 

1 



2 



average column density Navg ■ The median column density, 
or the column density most often be sampled, is 

Nme^ian = ^"^"^-'^ {N-Jj^-'^ + A^.-.g"^ V^/'""^^ (3) 

and in the limit of Nhigh = oo, N^edian = 2^/(™-^)7V;o^. 
Then, the ratio of the median to the mean is, in the limit 

of Nhigh = 00, 
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If the distribution function were infinitely sharply 
peaked, Nj^edian = -^2icm) but if the distribution has 
some significant width, these two quantities can be mea- 
surably different. For example, with a value of m = 3.5, 
consistent with the data of Wakker, Oosterloo & Putman 
(2002), N. median I N2icm = 0.792, Or about a 20% difference 
between the median and the mean. This implies that indi- 
vidual lines of sight out of the Galaxy will generally return 
a value for the column that is about 20% lower than the 
column induced from the 21 cm measurements. As an ex- 
ample of the importance of this in an astrophysical setting, 
this phenomenon will affect the interpretation of absorp- 
tion line measurements (single lines of sight) when they 
are compared to the amount of 21 cm emission (as was 
first pointed out to us by Wakker). In a comparison of a 
Lya absorption line to a 21 cm emission line measurement, 
one would find the systematic effect that the absorption 
column was less than the emission column. In the deter- 
mination of mctallicitics that compare metal absorption 
lines to 21 cm emission, one would determine a lower than 
true mctallicity more often than a higher metallicity. 

These small scale variations in the HI column will affect 
the absorption of light in the X-ray waveband in significant 
ways, especially at high optical depths. The transmitted 
fraction of light at some energy is given by 



Tr{E) = Io{E) 



J dx dy 



(5) 



At low optical depth {T{x,y,E) = a{E)N{x,y) « 1, 
where (j{E) is the energy dependent cross section), this 
just reduces to Tr{E) = Io{E){l - a{E)Navg), so it is ap- 
propriate to use Navg = N2icm, as is usually done. How- 
ever, if there are variations in N{x, y) when r > 1, the ef- 
fective X-ray absorption column changes with energy, with 
the deviation from A^2icm becoming greater toward lower 
energies, which is the direction of increasing optical depth. 
This can be illustrated in Figure 1, where we show the 
transmitted fractions as a function of energy for the case 



where Navq = N21C 



2 X 10^*^ cm ^ , and where the vari- 



ation occurs from Niow ~ 0.82 Navg to Nhigh = 1-31 Navg, 

with m = 3.5. At high optical depths, the transmission is 
underestimated when using A^2icm- The reason is due to 
the non-linearity in the behavior of the transmission with 
column, so that the fluctuations toward low density (lower 
optical depth) permit far more photons to flow through 
than an equivalent high density fluctuation removes. This 
leads to a residual flux when fltting a model with con- 
stant N2icm to the data, and this positive residual is a 
systematic effect that is prominent in a moderately nar- 
row energy range where 4 > r > 1. The absolute value 
of the residual is small at low energies because the opti- 
cal depth through the Galaxy is large, while the residual 



vanishes at larger energies because the optical depth be- 
comes low enough that using a constant N^icm is a good 
approximation. This effect can be rather important in the 
discussion of the reality of the soft excess X-ray emission, 
as we discuss below. However, first it is important to un- 
derstand the magnitude of the fluctuations in N so that 
the importance of the effects can be assessed. 

3. THE MAGNITUDE OF SMALL-SCALE STRUCTURE IN 
THE ISM FROM 21 CM STUDIES 

Most extragalactic objects that have been studied lie at 
Galactic latitudes b > 20° , or lines of sight that are out of 
the plane of the Milky Way. For such lines of sight, there 
are few studies of the small-scale structure in the total HI 
column, and almost no information on the structure of H2 . 
The type of information needed is the distribution func- 
tion of column densities within the larger beam of single- 
dish surveys (i.e., a histogram of N/f/). In principle, this 
can be obtained with synthesis arrays, but such instru- 
ments normally resolve away the larger scale structure, 
showing only small-scale variations. However, synthesis 
array observations, used in conjunction with single-dish 
observations (the zero spacing data) can recover all the 
emission and show the true distribution of N// / . Unfortu- 
nately, few such observations have been made out of the 
plane, although there are many such observations in the 
plane, which show a rich amount of structure. Out of the 
plane, two published observations are by Wakker, Oost- 
erloo & Putman (2002), which is toward a high velocity 
cloud (WW 187) in the direction of NGC 3783, / = 287°, 
h = 23°, and by Joncas, Boulanger & Dewdney (1992), 
which is toward a filament of far infrared cirrus emission 
at I = 141° b = 30°. Although these may not be typical 
regions, we examine them since they provide a starting 
point for understanding small-scale structure; we will be 
obtaining 21 cm data toward more typical sight lines in 
the near future. 

The observation by Joncas, Boulanger & Dewdney 
(1992) is toward a gas and dust complex in Ursa Major 
that has conspicuous IR-cirrus emission and is also known 
to have some molecular gas. However, the mean HI col- 
umn in this direction is 3.4x10^" cm^^ , which is fairly 
typical for sight lines at this Galactic latitude. The ob- 
servation was made with the Synthesis Radio Telescope 
of Dominion Radio Astrophysical Observatory, including 
zero spacing data and complementary observations in the 
continuum; the resulting resolution was 1' and the field 
of view was 2.6°. Since the sensitivity is dropping toward 
the edges, and the absolute calibration suffers, we have 
chosen to look at the variations in the central 1° region, 
smoothed to 2.3' and formed a histogram of the column 
densities (Figure 2). There is a fairly sharp cutoff in N(HI) 
for values below 1.7x10^'^ cm~^ , but there is a distribu- 
tion that extends to 5.9 x 10^*^ cm~^ , a range of almost a 
factor of four (0.54 in the log). 

The observation by Wakker is of a HVC projected upon 
a background AGN (NGC 3783) in the Southern Hemi- 
sphere, for the purpose of understanding the optical and 
ultraviolet absorption line results toward NGC 3783. The 
synthesis observations were made with the Australia Tele- 
scope Compact Array, supplemented with data from the 
Parkes Multibeam survey, so that HI structures were not 
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resolved out of the image. The map that they produced 
was not of all of the Galactic HI in that region, but only 
the HI in the velocity range surrounding cloud, or 190-280 
km s~^. Although it would be better to have the entire 
Galactic column, such a map was not produced, so wc have 
analyzed the existing map, which should give us some in- 
sight into HI variations. We have used a square box 1° on 
a side located in the central region in order to determine 
the column density histogram (Figure 3), which shows that 
there is a sharp drop in the cohmm below 0.65 x 10^*^ cm^^ 
. The HI distribution is moderately flat to 1.2 x 10^" cm~^ 
, declining toward higher values (16 x 10^° cm~^ in a way 
that may be described as a broken power-law of slopes -3 
to -4 (or with an exponential model). A detailed anal- 
ysis of this distribution is given in Wakker, Oosterloo & 
Putman (2002). 

Although there are few HI maps such as the above, the 
entire sky has been mapped in the far infrared with DIRBE 
and IRAS data by Schlegel, Finkbcincr & Davis (1998). 
These authors have produced maps of E(B-V), which are, 
on average, proportional to the HI column density, and 
with 2.67' resolution. We sampled five randomly chosen 
1° square regions located at l,b ^ -141.3°, 37.3°; 14.0°, 
74.8°; 112.8°, 47.7°; -44.3°45.4°; and 52.4°, 55.9°, and 
analyzed the column density distribution, using the rela- 
tionship of N{HI) = 8.18 X 10^^E{B -V) + 2.54 x 10^^ 
cm~^ , which was determined by us using other lines of 
sight out of the Galaxy. It is eviclcnit that the full range 
is at least 100% relative to the minimum value (Figure 4). 
However, we can better characterize the range by using the 
quartile points, referenced to the median column density, 
which in percentage, is 100 {N{75%) - N{25%))/Nmedian, 
where 7V(75%) and 7V(25%) arc the columns at the 75% 
and 25% quartile locations. Using this measure, we find 
that for these selected regions, the range is fairly narrow, 
14.6-19.2% with an average value of 16.5%. 

Unfortunately, it is nearly impossible to determine 
whether the degree of fluctuations seen in the HI is the 
same as in the IRAS-DIRBE data at a resolution of 1'. The 
HI data from Wakker show a greater degree of HI fluctu- 
ations than docs the IRAS-DIRBE data, but this may be 
because the HI map does not include HI at all velocities. 
In the region studied by Joncas, the HI and IRAS-DIRBE 
variations are similar in magnitude, but this region was 
chosen because of the prominent IRAS structure in the 
field, which guarantees substantial variation. One might 
expect the IRAS-DIRBE variations to be smaller than the 
HI variations just because the "beam" of the IRAS-DIRBE 
data is an order of magnitude larger than the HI radio 
synthesis beam of about 1' (so small scale variations are 
washed out). This is an issue where more data would be 
extremely valuable. 

There is additional evidence for small-scale structure in 
the interstellar medium, such as from the studies that 
are based upon optical and ultraviolet absorption mea- 
surements. The studies of Nal absorption systems toward 
globular cluster stars show that there are variations in the 
equivalent width by a factor of a few on scales of 1-4' 
(Langer, Prosser & Snedcn 1990) as well as on scales of 
~ 10" (Andrews, Meyer & Lauroesch 2001). Absorption 
line studies toward binary stars show that there is sig- 
nificant structure in the ISM on a scale of 0.1" (Laroesch, 



Meyer, and Blades 1999). In addition, HI absorption stud- 
ies of the 21 cm line show that 10-15% of cold HI occurs 
on scales below 0.1". Frail et al. (1994) study the changes 
in 21 cm absorption toward pulsars with known velocities, 
detecting changes in the absorption line over a timoscalo 
of a few years, which corresponds to a change of 10-100 
mas. On comparable angular scales, Faison et al. (1998) 
used VLBI studies toward radio AGNs and find variations 
in the HI optical depth as large as a factor of two. Unfor- 
timatcly, neither the optical absorption lines (Nal, KI) or 
the 21 cm absorption line is a true linear indicator of the 
gas column because the optical metal lines can be sensitive 
to changes in the ionization state while the 21 cm optical 
depth is inversely proportional to the temperature. Nev- 
ertheless, these studies are consistent with the presence 
of structure in the interstellar column on relatively small 
scales. 

4. SMALL-SCALE VARIATION FROM X-RAY OBSERVATIONS 

The absorption of soft X-rays depends upon the column 
of intervening gas as well as the energy of the observation, 
since the approximate relationship for the optical depth 
is r oc NE^'\ where is the column of X-ray absorbing 
gas (usually expressed as an equivalent HI column) and 
E is the energy of the photons. For a column density of 
TV = 2 X 10^° cm-2 , T = 1 at 0.24 keV, while for a col- 
umn of 4 X 10^" cm~^ , r = 1 occurs at 0.32 keV, so if 
there are significant fluctuations in the X-ray absorption 
column, we should see a "drop-out" of the soft emission 
(the ratio of soft to hard emission decreases). At these col- 
umn densities, the absorption of X-rays is due primarily to 
Hcl, with HI being the second most important absorber. 
Other components of the ISM are less important as ab- 
sorbers (warm ionized gas or hot X-ray emitting gas), so 
the X-ray column is a fairly good measure of the neutral 
column. 

This effect is, in principle, possible to detect against a 
relatively smooth extended extragalactic X-ray emitting 
object such as a cluster of galaxies. In practice, it is 
possible to search for this effect only in clusters of high 
surface brightness and for the correct range of Galactic 
column density, at least when using ROSAT data. The 
ROSAT PSPC data has rather limited energy resolution 
so that the data can only be separated into a two or three 
independent spectral regions. One of these divisions is due 
to the presence of carbon in the optical path, which leads 
to a sharp absorption edge at 0.284 keV that becomes pro- 
gressively less important at higher energies. Operationally, 
this leads to a local minimum in the detection efficiency 
near 0.5 keV, so wc have divided the PSPC data into two 
bands, 0.2-0.5 keV and 0.5-2.0 keV (PI channels 20-51 and 
52-201). Snowden et al. (1994) have developed procedures 
to extract diffuse emission in various energy bands, cor- 
rected for point sources, instrumental background, and lo- 
cal backgroimd sources (solar scattered X-rays, etc.). Here 
we use those procedures to form the ratio of band R2-I-R3 
(the 0.2-0.5 keV band) to R4-FR5-FR6-FR7 (0.5-2.0 keV). 

In order to detect variations in the ratio of these two 
bands, there must be sufficient photons in the soft band 
for statistical reasons and that implies a mean Galactic 
absorption column of less than about 6 x 10^° cm~^ (oth- 
erwise there are no soft photons) but the absorption in the 
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soft band should be significant (t = 1) so tliat modest vari- 
ations in r lead to detectable variations in the X-ray ratio. 
For example, if the mean Galactic column is 3 x 10^° cm~-^ 
and there is a region where the column rises to 4 x 10^*^ 
cm^^ , our hardness ratio will decrease by 30% while if 
there is a region where the column rises to 5 x 10^" cm~^ 
, the hardness ratio will drop to half its original value. 
Prom the above discussion, variations of this sort are not 
likely, although the probabilities of finding such variations 
can depend strongly on region. For the HI data of Joncas, 
Boulanger & Dewdney (1992) and Wakker, Oosterloo & 
Putman (2002) about 20% of the data are above the mean 
by 4/3 and 10% of the data are above the mean by a factor 
of 5/3 or greater. However, if one uses the variations in 
the DIRBE/IRAS data for these two regions and the five 
other regions (above), typically only 1% of the area has an 
intensity value that is above the mean by 4/3 or greater. 
Since the variation in the DIRBE/IRAS data is less than 
that in the HI for the region observed by Wakker, Oost- 
erloo & Putman (2002), this latter value of 1% may be 
an underestimate. This suggests that finding significant 
variation in our hardness ratio is likely if there are several 
hundred independent positions. 

We analyzed five clusters as given in Table 1, all of which 
are bright nearby systems that are fairly symmetrical in 
their inner 6'. The X-ray temperatures and the bolometric 
luminosities are taken from Wu, Xue & Fang (1999), for 
Hq = 75 km s^^ Mpc^^, and the HI columns are those 
given by Dickey & Lockman (1990). 

[EDITOR: INSERT TABLE 1 HERE.] 

Because galaxy clusters have radial temperatiirc gradi- 
ents, which would lead to radial softness variation, we ex- 
amined softness variations azimuthally, cither in circular 
annuli, or in the case of clusters with elliptical isophotes, 
in elliptical annuli. We avoided the central region and 
used four annular regions with radii (or semi-major axis) 
of l'-2' (annulus 1), 2'-3' (annulus 2), 3'-4.25' (annuhis 3) 
and 4.25'-5.5' (annulus 4). We also examined larger an- 
nuli, but the surface brightness was generally too low to 
be of use. In addition, each annular region was divided 
into 24 regions of 15° each, with region 1 being from po- 
sition angle 0° to 15°, with the position angle increasing 
east from north as is the usual convention. This leads to 
96 independent positions per cluster. 

The signature that we are looking for in the X-ray emis- 
sion is a local hardening of the emission due to the soft 
photons being absorbed by an excess in the Galactic HI. 
This will lead to a decrease in the hardness ratio with 
the change due to decrease in the number of soft photons 
rather than an increase in the number of hard photons. 
We discuss each of the systems in turn. 

Abell 119: The hardness ratio is shown as a function 
of annular ring and region (position angle) in Figure 5. 
There are three positions for which the hardness ratio is 
more than 3(t from the mean of the sample (annulus 1, 
region 2; annulus 2, region 19; annulus 3, region 8). In 
each case, the deviation is in the sense that the region is 
harder than the mean and it is the soft band that has de- 
creased relative to it neighbors rather than an increase in 
the hard band. These are precisely the expectations for 
regions with excesses in their Galactic HI columns. Also, 
the next four deviant points in the sample (about 2.5a; 



regions 6, 11, and 22 in annulus 3 and region 4 in annulus 
4) arc also harder than the average and due to a decrease 
in the number of soft photons. 

Abell 2029: There is one region that deviates by more 
than 3(7 and is in annulus 1 at region 10, but this is due to 
a positive deviation in the hard fiux rather than a decrease 
in the soft flux, so this does not meet our criteria. 

Abell 2142: In this system, we used elliptical isophotes 
with an eccentricity of 0.82 and an orientation of the major 
axis of —58°. There is a strong point source that coidd not 
be removed and it affects the first three regions in annuli 
2, 3, and 4 (at PA 0-45°), so they were removed from the 
analysis. The photon statistics were good for this long ob- 
servation, so the annuli were divided into 30 regions, each 
12° wide. There are two regions with deviations more 
than 3cr (annulus 2, region 26; annulus 3, region 8), both 
regions being harder than the average and the hardness is 
due to the lack of soft photons rather than an excess of 
hard photons (Figure 6). 

For the remaining two systems, AWM4 and Abell 3376, 
there was no evidence of deviations in the hardness ra- 
tio beyond 3a, although the numbers of photons were less 
than the other clusters, so the statistical significance is 
poorer. 

To summarize, there were about 600 locations in which 
we could have seen fluctuations in the hardness ratio, al- 
though about half those were of lower statistical quality 
than the others. Six loc;ations deviated from the mean 
hardness ratio in their annuli by more than 3a, all toward 
the region of fewer soft photons. About one such deviation 
is expected from statistical fluctuations, so this represents 
an excess. Five of these six regions were harder because 
of a paucity of soft photons, as would be expected if there 
were a positive excess in the Galactic neutral hydrogen 
column. We cannot rule out that these statistically mean- 
ingful fluctuations in the hardness ratio is intrinsic to the 
cluster. To determine whether the fluctuations are due to 
variations in the Galactic gas will require an accurate HI 
map in this region. In addition, we will be able to more ac- 
curately define the hardness variations with future X-ray 
observations, especially with XMM-Newton. 

Although point sources are largely removed in the data 
processing, we comment on the effect by point sources that 
are unresolved and therefore would not be removed. In the 
ROSAT band, the typical point source is of comparable 
hardness to the thermal Bremsstrahlung X-ray emission 
from a galaxy cluster. Also, the surface brightness of a 
galaxy cluster is much greater than that of the unresolved 
X-ray background sources. Taken together, we expect that 
unresolved point sources will not influence our analysis. 

5. APPLICATIONS 

The fluctuations in the foreground absorption column 
can have a number of effects, each of which should be 
modeled. These fluctuations can lead to variable absorp- 
tion and reddening, which will affect the properties of an 
HR diagram of a star cluster. Another example occurs 
in absorption line studies (as discussed above), since the 
mean neutral gas column inferred from 21 cm emission 
studies is likely to be higher than the median line of sight, 
which is more relevant for individual lines of sight. A third 
example, which we explore in more detail, is the modeling 
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of galaxy clusters, since this will be a systematic effect. 
This systematic effect either leads one to obtain an X-ray 
absorption column lower than the mean column (which 
would help explain the results of Arabadjis & Bregman 
(1999b), or alternatively, it will lead to a flux residual 
that might be attributed to a new "soft" component of 
the X-ray emission. 

The presence of a new soft X-ray component was first 
presented by Lieu et al. (1996) from data obtained with 
the EUVE satcUitc, combined with ROSAT X-ray data. 
The instrument used by the EUVE for this detection was 
the DS telescope (Boyer & Malina 1991), which is sensitive 
to photons up to about 0.19 keV, which is the soft X-ray 
region. Soft X-rays can be absorbed by Galactic gas, and 
for a moderately low column density of 2 x 10^" cm^^ , the 
optical depth to continuous absorption is unity at 0.24 keV 
and T = 3 at 0.16 keV, so the X-rays that are detected in 
the EUVE band lie in the heavily absorbed region of the 
X-ray spectrum. Lieu et al. (1996) took the X-ray spec- 
trum of some bright galaxy clusters and fit an absorption 
spectrum with a single column density and then compared 
the flux in the EUVE band to the measured value. They 
found that the EUVE emission was brighter than antici- 
pated, leading them to conclude that there was a very soft 
X-ray excess in the emission that needed to be explained. 

The validity of this emission has been questioned from 
two other works on the subject. Arabadjis & Bregman 
(2000) pointed out that the same soft excess should be 
visible in the ROSAT PSPC as well, and possibly easier 
to detect. The effective areas of the PSPC and the EUVE 
DS telescope are shown in Figure 7, which includes not 
only the raw collecting area of both instruments, but the 
effective; collecting are due to absorption by Galactic gas 
with a column density of 2 x 10^° cm~^ , which is fairly 
typical of the sight lines involved. One finds that the effec- 
tive EUVE bandpass is 0.144-0.186 kcV (FWHM), which 
can be a great aid because although the PSPC has more 
effective collecting area in this energy range, its energy res- 
olution is significantly poorer than the EUVE bandpass of 
0.04 keV. In this EUVE bandpass, the optical depth is 
about 3, so if the Galactic column has been overestimated 
slightly, it makes a large difference in the predicted flux. 
Galactic columns out of the plane, such as these, are dif- 
ficult to determine accurately, and Arabadjis & Bregman 
(1999b) showed that in nearly all cases, the soft X-ray ex- 
cess would vanish if one adopted a column that was less 
than 3cr from the published values. The only exception 
was the Coma cluster, for which the Galactic column is 
particularly small. For low HI columns, the corrections to 
the 21 cm observations are the greatest, so it is possible 
that the HI column is particularly low and that there is 
no soft X-ray excess. 

The validity of the soft X-ray excess was also questioned 
by Bowyer, Bcrghofer & Korpela (1999), who showed that 
certain data processing issues could complicate the abil- 
ity to analyze the data. In particular, they showed that 
flat-fielding was a key issue, so they constructed flat fields 
from many long observations, and when this was applied, 
soft excesses vanished in Abcll 2199 and Abell 1795, two 
clusters that Lieu, Bonamente & Mittaz (1999) claimed to 
have excesses. However, they find that the Coma cluster 
still has an excess in its soft X-ray emission, the same re- 



sult as Arabadjis & Bregman (1999b), working from the 
ROSAT data. 

Here we show that even if one has measured the mean 
Galactic 21 cm column correctly, fluctuations in that HI 
column lead to a soft X-ray excess. As an example, when 
we use the model described in §2, the residual flux has 
a peak at about 0.19 keV and a substantial amount of 
flux falls in the EUVE band (Figure 8). In the EUVE 
band, the fractional excess is claimed to be about 20-50% 
above the expected values, and that is approximately the 
amount that we predict from our model of column density 
fluctuations (Figure 9). 

Ideally, one should use the detailed column density map 
associated with the particular clusters in question, but 21 
cm data of this sort does not exist. Consequently, we 
cannot demonstrate that this effect is responsible for the 
extraction of a soft X-ray excess; we merely suggest it as 
a possibility. 

6. DISCUSSION AND FINAL COMMENTS 

We have shown the potential implications of variations 
in the Galactic absorption column in trying to interpret a 
variety of observations. The greatest shortcoming in deter- 
mining the true importance of this effect is the lack of accu- 
rate data on the HI variations on the sky at small angular 
scale. Such observations need a combination of synthesis 
array observations in combination with single dish data in 
order to include the "zero spacing" data. Aside from the 
two observations discussed above, data of this kind have 
been largely confined to the Galactic plane, where the Do- 
minion Radio Astronomy Observatory has undertaken a 
survey from 74° < I < 147°, but in the narrow latitude 
band —3.5° < b < 5.5°. This is far from the region where 
most of the galaxy clusters lie, typically 20° < 

Another issue that needs to be examined is the de- 
gree to which the HI fluctuations in these flelds follows 
the DIRBE/IRAS maps of Schlegel, Finkbeiner & Davis 
(1998). If it could be shown that the HI generally follows 
the DIRBE/IRAS data to good accuracy, tlicni one could 
use that data as a proxy for the HI data. Currently, the sit- 
uation is unclear since in the WW 187 field, there is a poor 
correlation between the HI column and the DIRBE /IRAS 
data, while in the IRAS filament mapped by Joncas, the 
correlation is much better. 

The are a variety of effects that will be clarified by fu- 
ture HI observations, so we strongly encourage observers 
to obtain such data. 

We are especially indebted to Bart Wakker and Gilles 
Joncas for sending us copies of their data for us to analyze. 
Also, we would like to thank Thomas Berghofer for his 
patience in answering our many questions about EUVE 
data. We would like to thank Jay Lockman, Jimmy Ir- 
win, Renato Dupke, Morton Roberts, and Stu Bowyer for 
their advice and comments. Partial support for this work 
has been provided by NASA through LTSA grant NAG5- 
10765. 
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Fig. 1. — The transmission fraction due to neutral atomic gas for a column of 2 X 10^'^ cm^'^ (dashed line) compared to the transmission for 
gas of the same mean column but with fluctuations as described in the text. The amount of light transmitted through the absorbing medium 
is less for a uniform absorption column then one that has fluctuations, so the failure to include such fluctuations will lead to a systematic 
residual at energies where the optical depth exceeds unity. 
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Fig. 2. — Histogram of the HI column density in the central 1° field that was obtained by Joncas et al. (1992). The few very low column 
density pixels is due to the S/N of the observation (an rms of 5 K). The range in the column density is more than a fa<;tor of three in this 
region, which includes a filament identified in IRAS data. 
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N(HI) (id"' cm"') 



Fig. 3. — Histogram of the HI column density in the field around WW 187 obtained by WaJjker, Oosterloo & Putman (2002). A power-law 
of slope -3.5 would be steeper than the observed distribution for 0.8 x lO^^cm"^ < N2icm. < 2 x lO^^cm"^, and it would be shallower than 
the observed distribution at greater columns. 
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Fig. 4. — The column density distribution within 1° fields inferred from the IRAS-DIRBE data of Schlegel, Finkbeiner & Davis (1998); 
arbitrary vertical shifts are introduced to separate the histograms. In general, there is a peak in the distribution, with a greater extension to 
higher column than to lower columns, and in some cases, there is a sharp cutofi' at the low column density side. 
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Fig. 5. — The ratio of the fluxes in the 0.2-0.5 keV band to that in the 0.5-2.0 keV band for 15° regions within annuli at l'-2' (annulus 1), 
2'-3' (annulus 2), 3'-4.25' (annulus 3), and 4.25'-5.5' (annulus 4), centered on the X-ray center of Abell 119. Several regions have hardness 
ratios that deviate signiflcantly from the mean, such as at: annulus 1, region 2; annulus 2, region 19; and annulus 3, region 8. Deviations of 
this nature can be caused by excess gas along the line of sight, which preferentially absorbs the soft X-ray emission. 
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Fig. 6. — The same ratio as in Figure 4, but for elliptical annuli centered on the X-ray center of Abell 2142. There is a strong source of 
extended emission going from annuli 2-4 in regions 1-3, so these regions were not used in calculating the mean hardness ratio. There are two 
regions that deviate more than 3a from the mean (annulus 2, region 26; annulus 3, region 8) and they are harder, consistent with excess 
absorption by foreground gas. 
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Fig. 7. — The raw and effective collecting areas for the EUVE DS telescope and the ROSAT PSPC in the energy range 0.11 - 0.23 IceV. 
The effective area is the product of the raw collecting area and the fractional transmission through a Galactic gas column of 2 X 10^" cm~'^ 
. The EUVE effective area has a bandwidth of about 0.14-0.19 keV (FWHM), so it samples a well-defined waveband, but it has an order of 
magnitude less effective collecting area than the ROSAT PSPC. 
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Fig. 8. — The residual between the incident flux and the model for the flux when using the mean column density (this is the diflference 
between the two lines in Figure 1). The dashed box shows the FWHM of the EUVE DS instrument along the energy axis. 
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Fig. 9. — The same residual shown in Figure 7, but expressed as tiic ratio of the residual flux to the model flux when using the average 
column density. Within the dashed box, which shows the effective bandwidth of the EUVE DS in energy, the residual flux can be 20-60% of 
the expected flux (a flux- weighted value of about 35%), comparable to the excess soft X-ray emission that has been discussed in the literature. 
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Table 1 
Galaxy Cluster Properties 



Cluster 


RA 


DEC 


z 


Lx 


T 


Nhi 


^cxp 


Name 


(J2000) 


(J2000) 




(10^4 erg s-i) 


(keV) 


(1020 cm-2 ) 


(ksec) 


Abell 119 


00:56:16.8 


-01:15:00 


0.044 


3.2 


5.6 


3.45 


15.2 


Abell 3376 


06:00:43.6 


-40:03:00 


0.0455 


2.1 


4.0 


4.81 


12.0 


Abell 2029 


15:10:58.7 


+05:45:42 


0.0767 


19. 


8.5 


3.18 


15.6 


Abell 2142 


15:58:16.1 


+27:13:29 


0.0899 


27. 


9.7 


3.86 


44.2 


AWM 4 


16:04:.")."). S 


+2;-!:55:.")4 


i)m2?> 


O.:-!!) 


2.4 


5.0;-! 


19.9 



